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Abstract —In this paper, we present a novel approach for robust 
optimal resource allocation with joint carrier aggregation to 
allocate multiple carriers resources optimally among users with 
elastic and inelastic traffic in cellular networks. We use utility 
proportional fairness allocation policy, where the fairness among 
users is in utility percentage of the application running on the 
user equipment (UE). Each UE is assigned an application utility 
function based on the type of its application. Our objective is 
to allocate multiple carriers resources optimally among users 
subscribing for mobile services. In addition, each user is guar¬ 
anteed a minimum quality of service (QoS) that varies based 
on the user’s application type. We present a robust algorithm 
that solves the drawback in the algorithm presented in 03 by 
preventing the fluctuations in the resource allocation process, in 
the case of scarce resources, and allocates optimal rates for both 
high-traffic and low-traffic situations. Our distributed resource 
allocation algorithm allocates an optimal rate to each user from 
all carriers in its range while providing the minimum price 
for the allocated rate. In addition, we analyze the convergence 
of the algorithm with different network traffic densities and 
show that our algorithm provides traffic dependent pricing for 
network providers. Finally, we present simulation results for the 
performance of our resource allocation algorithm. 

Index Terms —Optimal Resource Allocation; Joint Carrier Ag¬ 
gregation; Utility Proportional Fairness; Elastic Traffic, Inelastic 
Traffic 

I. Introduction 

In recent years, the number of mobile subscribers and their 
traffic have increased rapidly. Mobile subscribers are currently 
running multiple applications, simultaneously, on their smart 
phones that require a higher bandwidth and make users so 
limited to the carrier resources. Network providers are now 
offering multiple services such as multimedia telephony and 
mobile-TV More spectrum is required to meet these 
demands 0. However, it is difficult to provide the required 
resources with a single frequency band due to the scarcity of 
the available radio spectrum. Therefore, aggregating different 
carriers’ frequency bands is needed to utilize the radio re¬ 
sources across multiple carriers and allow a scalable expansion 
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of the effective bandwidth delivered to the user terminal, 
leading to interband non-contiguous carrier aggregation 0. 

Carrier aggregation (CA) is one of the most distinct features 
of 4G systems including Long Term Evolution Advanced (LTE 
Advanced). Given the fact that LTE requires wide carrier 
bandwidths to utilize such as 10 and 20 MHz, CA needs 
to be taken into consideration when designing the system 
to overcome the spectrum scarcity challenges. With the CA 
being defined in J5j, two or more component carriers (CCs) 
of the same or different bandwidths can be aggregated to 
achieve wider transmission bandwidths between the evolve 
node B (eNodeB) and the UE. An overview of CA framework 
and cases is presented in 0. Many operators are willing 
to add the CA feature to their plans across a mixture of 
macro cells and small cells. This will provide capacity and 
performance benefits in areas where small cell coverage is 
available while enabling network operators to provide robust 
mobility management on their macro cell networks. 

Increasing the utilization of the existing spectrum can 
significantly improve network capacity, data rates and user 
experience. Some spectrum holders such as government users 
do not use their entire allocated spectrum in every part of 
their geographic boundaries most of the time. Therefore, 
the National Broadband Plan (NBP) and the findings of the 
President’s Council of Advisors on Science and Technology 
(PCAST) spectrum study have recommended making the 
under-utilized federal spectrum available for secondary use 0. 
Spectrum sharing enables wireless systems to harvest under¬ 
utilized swathes of spectrum, which would vastly increase the 
efficiency of spectrum usage. Making more spectrum available 
can provide significant gain in mobile broadband capacity 
only if those resources can be aggregated efficiently with the 
existing commercial mobile system resources. 

This non-contiguous carrier aggregation task is a challeng¬ 
ing. The challenges are both in hardware implementation and 
joint optimal resource allocation. Hardware implementation 
challenges are in the need for multiple oscillators, multiple 
RF chains, more powerful signal processing, and longer bat¬ 
tery life j8j. In order to allocate different carriers resources 
optimally among mobile users in their coverage areas, a 
distributed resource allocation algorithm between the UEs and 
the eNodeBs is needed. 

A multi-stage resource allocation (RA) with carrier aggre¬ 
gation algorithms are presented in (9j-[lT)- The algorithm in 
0 uses utility proportional fairness approach to allocate the 


primary and the secondary carriers resources optimally among 
mobile users in their coverage area. The primary carrier first 
allocates its resources optimally among users in its coverage 
area. The secondary carrier then starts allocating optimal rates 
to users in its coverage area based on the users applications 
and the rates allocated to them by the primary carrier. A RA 
with CA optimization problem is presented in iflOll to allocate 
resources from the LTE Advanced carrier and the MIMO radar 
carrier to each UE, in a LTE Advanced cell based on the 
application running on the UE. A price selective centralized 
RA with CA algorithm is presented in llTTII to allocate multiple 
carriers resources optimally among users while giving the user 
the ability to select one of the carriers to be its primary carrier 
and the others to be its secondary carriers. The UE’s decision 
is based on the carrier price per unit bandwidth. However, 
the multi-stage RA with CA algorithms presented in m-an 
guarantee optimal rate allocation but not optimal pricing. 

In this paper, we focus on solving the problem of utility 
proportional fairness optimal RA with joint CA for multi¬ 
carrier cellular networks. The RA with joint CA algorithm 
presented in 12 fails to converge for high-traffic situations 
due to the fluctuation in the RA process. In this paper, we 
present a robust algorithm that solves the drawbacks in 12 
and allocates multiple carriers resources optimally among UEs 
in their coverage area for both high-traffic and low-traffic 
situations. Additionally, our proposed distributed algorithm 
outperforms the multi-stage RA with CA algorithms presented 
in 0-02 as it guarantees that mobile users are assigned 
optimal (minimum) price for resources. We formulate the 
multi-carrier RA with CA optimization problem into a convex 
optimization framework. We use logarithmic and sigmoidal- 
like utility functions to represent delay-tolerant and real-time 
applications, respectively, running on the mobile users’ smart 
phones m. Our model supports both contiguous and non¬ 
contiguous carrier aggregation from one or more network 
providers. During the resource allocation process, our dis¬ 
tributed algorithm allocates optimal resources from one or 
more carriers to provide the lowest resource price for the 
mobile users. In addition, we use a utility proportional fairness 
approach that ensures non-zero resource allocation for all users 
and gives real-time applications priority over delay-tolerant 
applications due to the nature of their applications that require 
minimum encoding rates. 

A. Related Work 

There has been several works in the area of resource 
allocation optimization to utilize the scarce radio spectrum 
efficiently. The authors in da-in have used a strictly 
concave utility function to represent each user’s elastic traffic 
and proposed distributed algorithms at the sources and the 
links to interpret the congestion control of communication 
networks. Their work have only focussed on elastic traffic 
and did not consider real-time applications as it have non¬ 
concave utility functions as shown in Oil- The authors in 
da and da have argued that the utility function, which 
represents the user application performance, is the one that 
needs to be shared fairly rather than the bandwidth. In this 


paper, we consider using resource allocation to achieve a utility 
proportional fairness that maximizes the user satisfaction. If a 
bandwidth proportional fairness is applied through a max-min 
bandwidth allocation, users running delay-tolerant applications 
receive larger utilities than users running real-time applications 
as real-time applications require minimum encoding rates and 
their utilities are equal to zero if they do not receive their 
minimum encoding rates. 

The proportional fairness framework of Kelly introduced 
in d3 does not guarantee a minimum QoS for each user 
application. To overcome this issue, a resource allocation 
algorithm that uses utility proportional fairness policy is 
introduced in H2- We believe that this approach is more 
appropriate as it respects the inelastic behavior of real-time 
applications. The utility proportional fairness approach in 
da gives real-time applications priority over delay tolerant 
applications when allocating resources and guarantees that 
no user is allocated zero rate. In 1321, (20) and (Dll . the 
authors have presented optimal resource allocation algorithms 
to allocate single carrier resources optimally among mobile 
users. However, their algorithms do not support multi-carrier 
resource allocation. To incorporate the carrier aggregation 
feature, we have introduced a multi-stage resource allocation 
using carrier aggregation in 0. In l22l and l23l . we present 
resource allocation with users discrimination algorithms to 
allocate the eNodeB resources optimally among mobile users 
with elastic and inelastic traffic. In (24), the authors have pre¬ 
sented a radio resource block allocation optimization problem 
using a utility proportional fairness approach. The authors 
in |25) have presented an application-aware resource block 
scheduling approach for elastic and inelastic adaptive real-time 
traffic where users are assigned to resource blocks. 

On the other hand, resource allocation for single cell multi¬ 
carrier systems have been given extensive attention in recent 
years If26) - l28l . In If29l - l32l . the authors have represented 
this challenge in optimization problems. Their objective is to 
maximize the overall cell throughput with some constraints 
such as fairness and transmission power. However, transform¬ 
ing the problem into a utility maximization framework can 
achieve better users satisfaction rather than better system¬ 
centric throughput. Also, in practical systems, the challenge is 
to perform multi-carrier radio resource allocation for multiple 
cells. The authors in (33), (34) suggested using a distributed 
resource allocation rather than a centralized one to reduce 
the implementation complexity. In l35l . the authors propose 
a collaborative scheme in a multiple base stations (BSs) 
environment, where each user is served by the BS that has 
the best channel gain with that user. The authors in [36) 
have addressed the problem of spectrum resource allocation 
in carrier aggregation based LTE Advanced systems, with the 
consideration of UEs MIMO capability and the modulation 
and coding schemes (MCSs) selection. 

B. Our Contributions 

Our contributions in this paper are summarized as: 

• We consider the RA optimization problem with joint CA 
presented in 12 that uses utility proportional fairness 


2 


approach and solves for logarithmic and sigmoidal-like 
utility functions representing delay-tolerant and real-time 
applications, respectively. 

• We prove that the optimization problem is convex and 
therefore the global optimal solution is tractable. In ad¬ 
dition, we present a robust distributed resource allocation 
algorithm to solve the optimization problem and provide 
optimal rates in high-traffic and low-traffic situations. 

• Our proposed algorithm outperforms that presented in m 
by preventing the fluctuations in the RA process when the 
resources are scarce with respect to the number of users. 
It also outperforms the algorithms presented in ®-ini 
as it guarantees that mobile users receive optimal price 
for resources. 

• We present simulation results for the performance of our 
RA algorithm and compare it with the performance of 
the multi-stage RA algorithm presented in 19I- 1TTI . 

The remainder of this paper is organized as follows. Section 
|II]presents the problem formulation. Section|nI]proves that the 
global optimal solution exists and is tractable. In Section [Tv] 
we discuss the conversion of the primal optimization prob¬ 
lem into a dual problem. Section [V] presents our distributed 
resource allocation algorithm with joint carrier aggregation 
for the utility proportional fairness optimization problem. In 
Section [VTl we present convergence analysis for the allocation 
algorithm and a modification for robustness at peak-traffic 
hours. In section I VIII we discuss simulation setup, provide 
quantitative results along with discussion and compare the 
performance of the proposed algorithm with the one presented 
in ED-QU- Section IVIIII concludes the paper. 

II. Problem Formulation 

We consider LTE mobile system consisting of AT carriers 
eNodeBs with I\ cells and M UEs distributed in these cells. 
The rate allocated by the I th carrier eNodeB to i th UE is 
given by ru where l = {1,2,...,A'} and i = {1,2,...,M}. 
Each UE has its own utility function 17* (r u + r 2 i + + rxi) 

that corresponds to the type of traffic being handled by the 
i th UE. Our objective is to determine the optimal rates that 
the I th carrier eNodeB should allocate to the nearby UEs. We 
express the user satisfaction with its provided service using 
utility functions that represent the degree of satisfaction of the 
user function with the rate allocated by the cellular network 
no na oa. We assume the utility functions Ui(ru + r 2 i + 
...+rKi) to be a strictly concave or a sigmoidal-like functions. 
The utility functions have the following properties: 

• Ui{ 0) = 0 and Uiiru + r 2i + ... + r K i) is an increasing 
function of ru for l. 

• Ui(ru+r 2 i+...+rKi) is twice continuously differentiable 
in ru for all l. 

In our model, we use the normalized sigmoidal-like utility 
function, as in m, that can be expressed as 

Ui{ru+ru + :.+r K i) = Ci( -E?-—-ci*) (!) 

V 1 _|_ e -Oi(E( = i r u-bi) ) 

where c,; = * and di = 1+e 1 a , b . . So, it satisfies Uj(0) = 

0 and Ui{ oo) = 1. We use the normalized logarithmic utility 


function, as in HE), that can be expressed as 


Ui(ru + r 2i + ... +r Ki ) 


iog(i + h Y?Ll r li) 

log(l + k max) 


( 2 ) 


where r max is the required rate for the user to achieve 100% 
utility percentage and ki is the rate of increase of utility 
percentage with allocated rates. So, it satisfies Ui{ 0) = 0 and 
Ui (rmax ) = 1- We consider the utility proportional fairness 
objective function that is given by 

M 

max n Ui (r u + r 2 i + ... + r K i) (3) 

2=1 


where r = {ri, r 2 ,..., r M } and r* = {m, r 2i ,..., r K i}- 
The goal of this resource allocation objective function is to 
maximize the total system utility while ensuring proportional 
fairness between utilities (i.e., the product of the utilities of 
all UEs). This resource allocation objective function inherently 
guarantees: 

• non-zero resource allocation for all users. Therefore, the 
corresponding resource allocation optimization problem 
provides a minimum QoS for all users. 

« priority to users with real-time applications. Therefore, 
the corresponding resource allocation optimization prob¬ 
lem improves the overall QoS for LTE system. 

The basic formulation of the utility proportional fairness 
resource allocation problem is given by the following opti¬ 
mization problem: 

M 

max n Ui(?u +r 2i + ... + rKi) 

2=1 

M M 

subject to E ru < Ai,y^r 2 j < R 2 ,... 

2=1 2=1 ^ 

M 

... , r Ki < Rk , 

i=l 

ru > 0, l = 1,2,..., AT, * = 1,2,..., M 


where Ri is the total available rate at the I th carrier eNodeB. 

We prove in Section [TIT] that the solution of the optimization 
problem (|4} is the global optimal solution. 


III. The Global Optimal Solution 
In the optimization problem {4}, since the objective func¬ 
tion argmaxnE Ui(r u + r 2i + ... + rKi) is equivalent to 

argmax^rfl! log(/7 i (ri l + r 2 i + ...+rif i )), then optimization 
problem Q can be written as: 


max E log (Ui(ru + r 2i + ... + r K i)) 

2=1 

M M 

subject to E ru < Ri, ^2 r 2 i < f? 2 , 

2=1 2=1 

M 

> E rKi — ^ k ' 

i= 1 

ru >0, l = 1,2 i = 1,2,..., M. 


( 5 ) 
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Lemma m.l. The utility functions \og(Ui(ru + ...in 
the optimization problem (0 are strictly concave functions. 

Proof: In Section HU we assume that all the utility 
functions of the UEs are strictly concave or sigmoidal-like 
functions. 

In the strictly concave utility function case, recall the utility 
function properties in Section UTl the utility function is positive 
Ui(ru + ... + rxi ) > 0, increasing and twice differentiable 


dUi{r 1 j + ...+r K i) 


and 9 2 u i (r li +...+r Ki ) 


d\og(Ui(ru + ... + r K i)) _ 


dru 


Ui 


> 0 


and 


d 2 \og(Ui(ru + ... + r Kl )) ^t Ui (dnT 


drl 


U? 


we have 


0 < C 7 ;(y 

di < 


1 


h \ di) K 1 
-p e ~ a i{I2i=i rii-bi) 

1 1 


< 


Cj d{ 


Ci 


1 q_ e -ai(T.f =1 r li-bi) 

— > 1 > 

di 

0 < 1 - dj(l + e ~ ai ^=i r a-bi)) < 


1 p Cidi 


1 


1 + Cidi 


d 

-— log Ufrii + ... + r Ki ) = 
UT'li 


dr u 


log Ufru + ... + r K i ) = 


Theorem III.2. The optimization problem 0 is a convex 
optimization problem and there exists a unique tractable 
global optimal solution. 


> 0 


with respect to ru. Then, it follows that --yy—- 

-jyyz -- < 0. It follows that, the utility function 


Proof: It follows from Lemma llll.ll that for all UEs utility 
functions are strictly concave. Therefore, the optimization 
problem © is a convex optimization problem |39l . The 
optimization problem 0 is equivalent to optimization problem 
0, therefore it is a convex optimization problem. For a convex 
optimization problem, there exists a unique tractable global 
optimal solution 


\og(Ui(ru + r 2 i + ... +rxi)) in the optimization problem 0 
have 


( 6 ) 


<0. (7) 


Therefore, the strictly concave utility function Ufru + r->% + 
... + rxi) natural logarithm logfjfm + r 2 i + ... + r K i)) 
is also strictly concave. It follows that the natural logarithm 
of the logarithmic utility function in equation 0 is strictly 
concave. 

In the sigmoidal-like utility function case, the util¬ 
ity function of the normalized sigmoidal-like function is 
given by equation 0 as Ufru + r 2 i + ... + r K i) = 


U:( - , — a. { , t~,— b,- ) - d i )• For 0 < EfcLlHi < 


It follows that for 0 < Ed =i Hi < Ei=l we ^ ave tFe first 
and second derivative as 


1 — di(l + e~ a dT,f=i 

aie~ ai ^h= i 

-i-:- > n 

— a- die~ a dT.f^x ni-bi) 
d (l -di{ 1 + e- a (S£i ' 

- fl ? e -«i(Ef=inry 

-|- 1 - < 0 

(1 _p e -a,i(J2f=i ru-bi)y 


IV. The Dual Problem 

The key to a distributed and decentralized optimal solution 
of the primal problem in 0 is to convert it to the dual problem 
similar to m, in and |40) . The optimization problem 0 
can be divided into two simpler problems by using the dual 
problem. We define the Lagrangian 

M 

L( r, p) = ^2 log (Ui(m + r 2l + ... + r K %) 

2=1 

M 

-pi(^2 r u+z i - i?i) -... 

2=1 
M 

— Pk(^ 2 TKi + Zk — Rk) 

2=1 
M 

-- ^ ( lo g{Ui(ru + r 2 i + ... + r K i)) ~ ^Pini) 
2=1 1 = 1 

K 

5 Z p i( R i - Zl ^ 

1=1 

4 K 

= +^2Pl(Rl ~ z l) 

2=1 1 = 1 

( 8 ) 

where zi > 0 is the I th slack variable and pi is Lagrange 
multiplier or the shadow price of the I th carrier eNodeB (i.e. 
the total price per unit rate for all the users in the coverage area 
of the I th carrier eNodeB) and p = {pi,p 2 , ■ ■■■!> k}- Therefore, 
the i th UE bid for rate from the I th carrier eNodeB can be 
written as w u = Pirn and we have Ei=i w u = Pi Eiti r li- 
The first term in equation 0 is separable in r,. So we have 
m r axE^(log(EMni + r 2 i + ■■■ + r K i)) - E?=iPi r u) = 
EZi max(log(C/j(rii + r 2i + ... +r K i)) ~ Eti Pim)- 
The dual problem objective function can be written as 


K 


M 


D( p) =max L(r, p) 


M 


K 


= y^max(L i (r i ,p)) - zf) 


(9) 


i=i 


Therefore, the sigmoidal-like utility function Ufru+.-.+rxi) 
natural logarithm \og(Ui(ru + ... + rKi)) is strictly concave 
function. Therefore, all the utility functions in our model have 
strictly concave natural logarithm. ■ 


The dual problem is given by 

D( p) 


min 

p 


subject to pi > 0, l = 1,2,..., K. 


( 10 ) 
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So we have 


dD(p) 

dpi 


M 


= Ri-Y, ru ~ zi = 0 


(id 


substituting by YJiLi w U = Pi r U we have 


Tn =1 Wli 

Ri - zi 


pi = 


( 12 ) 


Now, we divide the primal problem © into two simpler 
optimization problems in the UEs and the eNodeBs. The i th 
UE optimization problem is given by: 


max 

f'i 

subject to 


K 

log(£/i(rii + r 2i + ... + r Ki )) - ^pirn 

1=1 (13) 

Pi> 0 

ru> 0, i = l,2,...,M,l = l,2,...,K. 


The second problem is the I th eNodeB optimization problem 
for rate proportional fairness that is given by: 


min D(p) 

pi 

subject to pi > 0. 


(14) 


The minimization of shadow price pi is achieved by the 
minimization of the slack variable zi > 0 from equation 
(fl2l >. Therefore, the maximum utility percentage of the I th 
eNodeB rate Ri is achieved by setting the slack variable 
zi = 0. In this case, we replace the inequality in primal 
problem © constraints by equality constraints and so we have 
w n = PiRi- Therefore, we have pi = w ‘ z where 
wu = Pifu is transmitted by the i th UE to I th eNodeB. The 
utility proportional fairness in the objective function of the 
optimization problem © is guaranteed in the solution of the 
optimization problems dT3l) and ( 11 4b . 


V. Distributed Optimization Algorithm 

The distributed resource allocation algorithm, in ID, for 
optimization problems <o and d is a modified version 
of the distributed algorithms in lfl2 . [f20j . |2P . Ifl3l and 
SOI, which is an iterative solution for allocating the network 
resources for a single carrier. The algorithm in ID allocates 
resources from multiple carriers simultaneously with utility 
proportional fairness policy. The algorithm is divided into the 
i th UE algorithm as shown in Algorithm 1 © and the I th 
eNodeB carrier algorithm as shown in Algorithm 2 ©. In Al¬ 
gorithm 1 and 2 ©, the i th UE starts with an initial bid wu(l) 
which is transmitted to the I th carrier eNodeB. The I th eNodeB 
calculates the difference between the received bid wuin) and 
the previously received bid wu(n — 1) and exits if it is less 
than a pre-specified threshold S. We set wu{ 0) = 0. If the 
value is greater than the threshold, the I th eNodeB calculates 

' M / \ 

the shadow price pi ( n ) = w and sends that value to 

all UEs in its coverage area. The i th UE receives the shadow 
prices pi from all in range carriers eNodeBs and compares 
them to find the first minimum shadow price Pmin( n ) ar, d th e 
corresponding carrier index l\ £ L where L = {1,2,..., K). 
The i th UE solves for the l\ carrier rate that maximizes 



Fig. 1. Flow Diagram with the assumption that the shadow price from the 
first carrier eNodeB pi is less before the nith iteration so rate rn of the i th 
user is allocated. After the nith iteration, the shadow price from the second 
carrier eNodeB p 2 is less so rate r^i is allocated. 


log Ui{ru + ... + r Ki ) - J2iLiPi( n ) r n with respect to r hi . 
The rate rj(n) = ri x i{n) is used to calculate the new bid 
wi x i(n) = Pmi n (n)rj(n). The i th UE sends the value of its 
new bid wi^n) to the l\ carrier eNodeB. Then, the i th UE 
selects the second minimum shadow price p^ lin (n) and the 
corresponding carrier index l 2 £ L. The i th UE solves for 
the l 2 carrier rate ri 2 i(n ) that maximizes \ogUi(ru + ... + 
rio) - Ya=\ Pi{n)ru with respect to n 2i . The rate n 2 i {n) 
subtracted by the rate from l\ carrier rf(n) = ri 2 i(n ) —rj(n) 
is used to calculate the new bid wi 2 i(n ) = p{ rlin (n)r^(n) 
which is sent to l 2 carrier eNodeB. In general, the i th UE 
selects the m th minimum shadow price p™ in (n) with carrier 
index l m £ L and solves for the l m carrier rate ri m i(n) that 
maximizes logUi(ru + ...+rKi) — '£2p=iPi( n ) r ii with respect 
to n mi . The rate r; mi (n) subtracted by i carriers 

rates r™(n) = n mi (n) - (rj(n) + rf(n) + ... + rj n_ 1 (n)) 
is used to calculate the new bid wi m i(n) = p™ in (n)r™ (n) 
which is sent to l m carrier eNodeB. This process is repeated 
until | wu{n) — wu(n — 1)| is less than the threshold S for all 
l carriers. 

The distributed algorithm in ID is set to avoid the situation 
of allocating zero rate to any user (i.e. no user is dropped). 
This is inherited from the utility proportional fairness policy 
in the optimization problem, similar to Q2, ED and ED- In 
addition, the UE chooses from the nearby carriers eNodeBs 
the one with the lowest shadow price and starts requesting 
bandwidth from that carrier eNodeB. If the allocated rate is 
not enough or the price of the resources increases due to high 
demand on that carrier eNodeB resources from other UEs, the 
UE switches to another nearby eNodeB carrier with a lower 
resource price to be allocated the rest of the required resources. 
This is done iteratively until an equilibrium between demand 
and supply of resources is achieved and the optimal rates are 
allocated in the LTE mobile network. Figure © shows a block 
diagram that represents the distributed RA algorithm. 
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VI. Convergence Analysis 

In this section, we present the convergence analysis of 
Algorithm 1 and 2 in JT) for different values of carriers 
eNodeBs rates Hi. This analysis is equivalent to low and 
high-traffic hours analysis in cellular systems (e.g. change in 
the number of active users M and their traffic in the cellular 
system lf20ll ). 


A. Drawback in Algorithm 1 and 2 in [T] 

Lemma VI. 1. For sigmoidal-like utility function Ui(ru+r 2 i + 
...+rpci), the slope curvature function 9 log u d r ^t+^+..-+rKt) 
has an inflection point at E/Li r ii = r® ~ bi and is convex 

M E/=i r H > rf. 


Proof: For the sigmoidal-like function Ui(ru 


tnction Ui(ru + r^i + 
- let Si(m) = 


■■■ + r K i) = Ci[ - _ a ,^K ri . _ b , 

V 1+e l lx 

d\ogU i (r li +r 2 i +...+r Ki ) ^ s j 0 p e curva t ure function. Then, 
we have that 


dS^ 

dru 


—a^die a dJ2i=irn N) 
QT? g — a i(^2l = l ru-bi) 


4 - e -ai(Ei=i r H -bi)y 


and 

d 2 S z 

K 


x 3rf. e -a i (E^ 1 n i -6 i )(l _ d .(l _ e-adZfLi ru-bi))) 

Ci(l - di{ 1 + e -“'Ef=i r H-W)J 3 
a ^e~ a dZiLi r ii- bi )^i — e - a i(Z!Li r n- b A) 

(l -|_ e -Oi(£i=i r,i-bi)\ 3 

(15) 

We analyze the curvature of the slope of the natural logarithm 
of sigmoidal-like utility function. For the first derivative, we 
have <0 V ru- The first term S'} of in equation 


m can be written as 

i _ 0 f i e aibi {e a ' bi + e-^^Dn-bi)) 

( e €Libi — g — = i r H — ^3 

and we have the following properties: 

( lim Zf =1 ru^0 S } = °°’ 


lim 


ZlUm^bi 


Sj = 0 for bi > -1-. 


(16) 


(17) 


For second term Sf of 2—if. i n equation ( IT5l) . we have the 



(18) 


Sf(ru = h — Yttfi fji) = 0, 

Sfim > bi - Yij^Li r ji) > 0, 

Sf(rn <b z - rji) < 0. 

From equation (fTTb and (IT8l) . S', has an inflection point at 


!la=i r n = r t ~ bi . In addition, we have the curvature of Si 
changes from ; 
function befon 

convex function after the inflection point. 


changes from a convex function close to origin to a concave 
function before the inflection point Ez=i r u = r t then to a 


Our rate allocation approach guarantees non-zero rate allo¬ 
cation for all active users in the coverage area of a specific 
carrier eNodeB. We define the set AA l := {i : ru ^ 0} to be 
the set of active users covered by the I th eNodeB. Then, we 
have the following Corollary. 

Corollary VI.2. If r T f < If V / e L then Algo¬ 

rithm 1 and 2 in SD converge to the global optimal rates 
which correspond to the steady state shadow price p ss < 

a.maxtkma.x _|_ £l>ga* W /j ere j max = arg maX^g^I bi. 

Proof: For the sigmoidal-like function Ui(ru + r 2 i + 
• + r Ki ) = c i( e -a ig 1 R_ iT ., i -i, i) - d i)’ the optimal so¬ 

lution is achieved by solving the optimization problem 
In Algorithm 1 Q], an important step to reach to the opti¬ 
mal solution is to solve the optimization problem ru(n) = 
argmax^log f/j(ru + r 2i + ... + r K i) - Pi{n)nfj for every 

UE in the I th eNodeB coverage area. The solution of this 
problem can be written, using Lagrange multipliers method, 
in the form 


(9 log Ujir-q + r 2i + ■■■ + r Ki ) 
dru 


-pi = Si(ru) ~Pi = 0. 


(19) 


From equation (ITTb and (fill) in Lemma 1 VI. 1 1 we have the 
curvature of Si(ru ) is convex for Yffi r u > rf « bi. The 
algorithm in £0 is guaranteed to converge to the global optimal 
solution when the slope S’j (ju ) of all the utility functions 
natural logarithm logi/* (ru + r 2 »+... + r^i ) are in the convex 
region of the functions, similar to analysis of logarithmic 
functions in IE) and l40l . Therefore, the natural logarithm of 
sigmoidal-like functions log Ui(ru+r 2 i + ■■■ + rKi ) converge 
to the global optimal solution for Ez*=i r u > r® ~ bi. 
The inflection point of sigmoidal-like function Ui(ru + r 2 i + 
... + r Ki ) is at rf = bi. For r f < R h the 

algorithm in m allocates rates Ez^i r n > bi for all users. 
Since Si(ru ) is convex for Ezli r u > rf « bi then the 
optimal solution can be achieved by Algorithm 1 and 2 in 
m. We have from equation ( IT9t and as Si(ru) is convex for 
E*i r n > r ! ~ bi , that p ss < Si(J2iL i ru = max l£A1 i bi) 
where r u = m &x ieM i bi) = a \ 1 £; dimax + and 

± “ l max 

(max = argmax, eJM i b t . ■ 

We define the set M c := {i : ru f 0 V l £ £, ru = 0 V l £ 
C} to be the set of active users covered exclusively by the 
set of carriers eNodeBs C C L. Then, we have the following 
Corollary. 

Corollary VI.3. For ExeA/L r f > Eze£ r i an d the global 
optimal shadow price p ss « — aidie 2 a . b - —|— aie a 2 . t . where 

l—di(iye~2~ L ) (1+e-V*-) 

i £ M c , then the solution given by Algorithm 1 and 2 in ]/ 
fluctuates about the global optimal rates. 


Proof: For the sigmoidal-like function Ui(ru 4- r 2 i + 
■ ■■ + r K i) = c z y—-—— - dfj, it follows from 

lemma |VI.1| that for YU^M C r 'f > Ez G £ ^ 3 • y \T £ 

such that the optimal rates ErLi r h Pt < bi. Therefore, if 

p ss « — aidie 2 a . h - . —|— a,e a ) b . is the optimal shadow price 

l-ddl+e-TF) (1+eT) 
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Algorithm VI.l The i th UE Algorithm 

Send initial bid wu( 1) to I th carrier eNodeB (where l £ 

L = {1,2, K}) 

loop 

Receive shadow prices pieiA' 1 ) from all in range carriers 
eNodeBs 

if STOP from all in range carriers eNodeBs then 
Calculate allocated rates r? pt = Wli f n ) 

h Pl (n) 

STOP 

else 

Set Pmi„ = {} and rf = 0 
for m = 1 —>• I\ do 

Pmin(^) m in(p \ {Pmin’Pmini ■••iPmin }) 
l m = {l G L : Pi = min(p \ 
{Pmin^min^-^min 1 })} { 1 m « the index of the 
corresponding carrier} 

Solve ri i(n) = argmaxf log Uflru + ... + r K i) - 
r l m i V 

"Yld=iPl{ n ) r li\ f° r the l rn carrier eNodeB 

r?(n) = ri mi {n) - E™}, 1 rj (n) 
if r™(n) < 0 then 
Set rf (n) = 0 

end if 

Calculate new bid wi m i(n ) = Pmin( n ) r T ( n ) 
if | wi m i(n) — u>i m i(n — 1)| > A w(n) then 

wi m i(n) = Wi(n— l) + sign (wi mi (n) - wi mi (n- 
l))Aw(n) (An; = h\e ^ or Aw = 

end if 

Send new bid wi m i{n) to l m carrier eNodeB 

end for 
end if 
end loop 


for optimization problem (0). Then, a small change in the 
shadow price pi{n) in the n th iteration can lead the rate ru(n) 
(root of Si(ru) — pi{n) = 0) to fluctuate between the concave 
and convex curvature of the slope curve Sflru ) for the i th 
user. Therefore, it causes fluctuation in the bid wu (n) sent to 
the eNodeB and fluctuation in the shadow price pi(n) set by 
eNodeB. Therefore, the iterative solution of Algorithm 1 and 
2 in III fluctuates about the global optimal rates i r °T- ® 

Theorem VI.4. Algorithm 1 and 2 in HU does not converge 
to the global optimal rates for all values of Ri. 

Proof: It follows from Corollary IVI.2I and IVI.3I that 
Algorithm 1 and 2 in 0] does not converge to the global 
optimal rates for all values of Ri. ■ 

B. Solution usins Alsorithm \VI.l\ and WI.2\ 

For a robust algorithm, we add a fluctuation decay function 
to the algorithm presented in Q]] as shown in Algorithm lVI.il 
Our robust algorithm ensures convergence for all values of 
the carriers eNodeBs maximum rate Ri for all l. Algorithm 
IVI.ll and IVI.2I allocated rates coincide with Algorithm 1 and 
2 in m for JAieM 1 r f < Ri V l G L. For r f > 


Algorithm VI.2 The I th eNodeB Algorithm 

loop 

Receive bids wu(n ) from UEs {Let wu(0) = 0 Vi} 
if \wn(n) — wn(n — 1)| < S Vi then 
Allocate rates, r° pt = to i th UE 

’ li pi(n) 

STOP 

else 

Calculate pi(n) = 

Send new shadow price pi (n) to all UEs 

end if 
end loop 


Y^irc Ri' robust algorithm avoids the fluctuation in the non- 
convergent region discussed in the previous section. This 
is achieved by adding a convergence measure A w(n) that 
senses the fluctuation in the bids In case of fluctuation, it 
decreases the step size between the current and the previous 
bid wu(n) — wu(n — 1) for every user i using fluctuation 
decay function. The fluctuation decay function could be in the 
following forms: 

« Exponential function : It takes the form Aw(n) = 

(iie _7r 2. 

• Rational function: It takes the form A w(n) = —. 
where can be adjusted to change the rate of decay 

of the bids wu. 

Remark VI.5. The fluctuation decay function can be included 
in the UE or the eNodeB Algorithm. 

In our model, we add the decay part to the UE Algorithm 
as shown in Algorithm lVI.il 

VII. Simulation Results 

Algorithm IVI.ll and IVI.2I were applied to various logarith¬ 
mic and sigmoidal-like utility functions with different parame¬ 
ters in MATLAB. The simulation results showed convergence 
to the global optimal rates. In this section, we present the 
simulation results for two carriers in a heterogeneous network 
(HetNet) that consists of one macro cell, one small cell and 
12 active UEs as shown in Figure [2] The UEs are divided into 
two groups. The 1 st group of UEs (index i = {1,2, 3,4, 5, 6}) 
is located in the macro cell under the coverage area of 
both the 1 st carrier (Cl) and the 2 nd carrier (C2) eNodeBs. 
We use three normalized sigmoidal-like functions that are 
expressed by equation ([I]) with different parameters. The used 
parameters are a = 5, b = 10 corresponding to a sigmoidal- 
like function that is an approximation to a step function at 
rate r = 10 (e.g. VoIP) and is the utility of UEs with indexes 
i = {1,7}, a = 3, b = 20 corresponding to a sigmoidal- 
like function that is an approximation of an adaptive real¬ 
time application with inflection point at rate r = 20 (e.g. 
standard definition video streaming) and is the utility of UEs 
with indexes i = {2,8}, and a = 1, b = 30 corresponding 
to a sigmoidal-like function that is also an approximation of 
an adaptive real-time application with inflection point at rate 
r = 30 (e.g. high definition video streaming) and is the utility 
of UEs with indexes i = {3,9}, as shown in Figure 0 We 


7 































App of 2 



Fig. 2. System model with two groups of users. The 1 st group with 
UE indexes i = {1,2,3,4,5,6}, 2 nd group with UE indexes i = 
{7,8,9,10.11,12}. 


Uiiru + r 2 t) 



Fig. 4. The sigmoidal-like utility Ui (ru +r 2 i ) = a ( 1 +e -a l ( rii +r- 2 i -b l ) - 

di) of the i th user, where ru is the rate allocated by 1 st carrier eNodeB 
and V‘ 2 i is the rate allocated by 2 nd carrier eNodeB. 


TABLE I 

Users and their applications utilities 



Fig. 3. The users utility functions Uiiru + r 2 f) used in the simulation 
(three sigmoidal-like functions and three logarithmic functions). 

use three logarithmic functions that are expressed by equation 
(0 with r max = 100 and different fc; parameters which are 
approximations for delay-tolerant applications (e.g. FTP). We 
use k = 15 for UEs with indexes * = {4, 10}, k = 3 for UEs 
with indexes i = {5,11}, and k = 0.5 for UEs with indexes 
i = {6,12}, as shown in Figure [3 A summary is shown in 
table HI A three dimensional view of the sigmoidal-like utility 
function Ui(ru + r 2 i) is show in Figure 0] 

A. Allocated Rates for 30 < R\ < 200 and R 2 = 70 

In the following simulations, we set 5 = 10 -3 , the 1 st 
carrier eNodeB rate R\ takes values between 30 and 200 with 
step of 10, and the 2 nd carrier eNodeB rate is fixed at R 2 = 70. 
In Figure [3 we show the final allocated optimal rates r, = 
r 1 i + r 2i of different users with different I s * carrier eNodeB 
total rate Ri and observe how the proposed rate allocation 
algorithm converges when the eNodeBs available resources are 


Applications Utilities Parameters 

Users Indexes 

Sigi 

Sig a = 5, b = 10 

< = {1.7} 

Sig2 

Sig a = 3, b = 20 

i = {2,8} 

Sig3 

Sig a = l, b = 30 

i = {3,9} 

Logl 

Log k = 15, r m ax = 100 

i = {4,10} 

Log2 

Log h = 3, Tmax — 100 

i = {5,11} 

Log3 

Log k = 0.5, rmax = 100 

i = {6,12} 


abundant or scarce. In Figure [5(a)] we show the rates allocated 
to the 1 st group of UEs by only Cl eNodeB since C2 eNodeB 
is not within these users range, we observe the increase in the 
rate allocated to these users with the increase in R\. Figure 
|5(b)| shows the final allocated rates to the 2 nd group of UEs 
by both Cl and C2 eNodeBs. Since these users located under 
the coverage area of both the macro cell and the small cell, 
they are allocated rates jointly using the proposed RA with 
joint CA approach. Figure |5(a)| and |5(b)| show that by using 
the RA with joint CA algorithm, no user is allocated zero 
rate (i.e. no user is dropped). However, the majority of the 
eNodeBs resources are allocated to the UEs running adaptive 
real-time applications until they reach their inflection rates 
the eNodeBs then allocate more resources to the UEs with 
delay-tolerant applications, as real-time application users bid 
higher than delay-tolerant application users by using the utility 
proportional fairness policy. 

In Figure [6] we show the rates allocated to the 2 nd group 
users, located under the coverage area of both the macro cell 
and small cell eNodeBs, by each of the two carriers’ eNodeBs 
with the increase in the 1 st carrier eNodeB resources. In Figure 
|6(a)| and |6(b)[ when the resources available at C2 eNodeB (i.e. 
R 2 ) is more than that at Cl eNodeB, we observe that most of 
the 2 nd group rates are allocated by C2 eNodeB. However, the 
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(a) The rates allocated ru from the 1 st carrier eNodeB (i.e. the macro cell 
eNodeB) to users of the 1 st group (i.e. i = 1,2,3,4, 5,6). 



Ri 


(a) The allocated rates ru from the 1 st carrier eNodeB to the 2 nd group of 
users. 




(b) The rates ru + r 2 i allocated from 1 st and 2 nd carriers eNodeBs (i.e. (b) The allocated rates r 2 i from the 2 nd carrier eNodeB to the 2 nd group of 

the macro cell and the small cell eNodeBs) to users of the 2 nd group (i.e. users. 

i = 7,8,9,10,11,12). 

Fig. 6. The allocated rates from Cl and C2 eNodeBs to the 2 nd group of 
Fig. 5. The allocated rates ru of the two groups of users verses 1 st users with l s ‘ carrier eNodeB rate 30 < Ri < 200 and 2 nd carrier eNodeB 

carrier rate 30 < Ri < 200 with 2 nd carrier rate fixed at R 2 = 70. rate fi xe d at R 2 = 70. 


delay tolerant applications are not allocated much resources 
since most of Ri is allocated to the real-time applications. 
With the increase in Cl eNodeB resources f?i, we observe a 
gradual increase in the 2 nd group rates allocated to real-time 
applications from Cl eNodeB and a gradual decrease from 
C2 eNodeB resources allocated to real-time-applications. This 
shift in the resource allocation increases the available resources 
in C2 eNodeB to be allocated to 2 nd group delay tolerant 
applications by C2 eNodeB. 

B. Pricing Analysis and Comparison for 30 < R\ < 200 and 
Ri = 70 

In the following simulations, we set 5 = 10 -3 and the 1 st 
carrier eNodeB rate R± takes values between 30 and 200 with 
step of 10, and C2 eNodeB total rate is fixed at Ri = 70. As 
discussed before, the users’ allocated rates are proportional 
to the users’ bids. Real-time application users bid higher than 
delay-tolerant application users due to their applications nature 
and the utility proportional fairness policy. Therefore, the 
pricing which is proportional to the bids is traffic-dependent, 
i.e. when the demand by users increases, as a result the price 
increases and vice versa. 


In Figure [7] we compare between the shadow price of Cl 
and C2 eNodeBs when using the proposed RA with joint CA 
approach with their shadow prices obtained when using the 
multi-stage RA with CA approach in |9l- lfTTI . For the RA 
with joint CA case, we observe that the shadow price of Cl 
eNodeB is higher than that of C2 eNodeB for R\ < 80 
and approximately equal for 80 < R\ < 200 which shows 
how it is very efficient to use the joint CA approach for the 
pricing of the user. We also show how the prices decrease 
with the increase in the eNodeBs total rate. By using this 
traffic-dependent pricing, the network providers can flatten the 
traffic specially during peak hours by setting traffic-dependent 
resource price, which gives an incentive for users to use the 
network during less traffic hours. On the other hand, for the 
multi-stage RA with CA approach, we show in Figure [7] the 
changes in Cl and C2 eNodeBs shadow prices with R \. When 
using the multi-stage RA with CA approach, all users are first 
allocated rates by the macro cell eNodeB, once Cl eNodeB is 
done allocating its resources C2 eNodeB starts allocating its 
resources only to the 2 nd group users as they are located within 
its coverage area. Since the pricing method in multi-stage 
RA with CA approach is not optimal, this explains why the 
shadow prices of Cl and C2 eNodeBs, in Figure [7] when using 
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Fig. 7. The 1 st carrier shadow price pi and 2 nd carrier shadow price p 2 
for both multi-stage RA with CA and joint RA methods with Cl eNodeB rate 
30 < Ri < 200 and C2 eNodeB rate Ri = 70. 


the proposed RA with joint CA approach are less than their 
corresponding prices when using the multi-stage RA with CA 
approach. This shows how the proposed algorithm outperforms 
the algorithms presented in m-m as it guarantees that 
mobile users receive optimal price (minimum) for resources. 

VIII. Conclusion 

In this paper, we introduced a novel resource allocation 
optimization problem with joint carrier aggregation in cellular 
networks. We considered mobile users running real-time and 
delay-tolerant applications with utility proportional fairness 
allocation policy. We proved that the global optimal solution 
exists and is tractable for mobile stations with logarithmic and 
sigmoidal-like utility functions. We presented a novel robust 
distributed algorithm for allocating resources from different 
carriers optimally among the mobile users. Our algorithm 
ensures fairness in the utility percentage achieved by the 
allocated resources for all users. Therefore, the algorithm gives 
priority to users with adaptive real-time applications while 
providing a minimum QoS for all users. In addition, the 
proposed RA with joint CA algorithm guarantees allocating 
resources from different carriers with the lowest resource price 
for the user. We analyzed the convergence of the algorithm 
with different network traffic densities and presented a robust 
algorithm that overcomes the fluctuation in allocation during 
peak traffic hours. We showed through simulations that our 
algorithm converges to the optimal resource allocation and 
that the proposed algorithm outperforms the multi-stage RA 
with CA algorithms presented in m-m as it guarantees that 
mobile users receive optimal price for the allocated resources. 
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